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Abstract 

We derive a path integral expression for the transition amplitude in 1+1-di- 
mensional QCD starting from canonically quantized QCD. Gauge fixing after quan- 
tization m leads to a formulation in terms of gauge invariant but curvilinear vari- 
ables. Remainders of the curved space are Jacobians, an effective potential, and sign 
factors just as for the problem of a particle in a box. Based on this result we de- 
rive a Faddeev-Popov like expression for the transition amplitude avoiding standard 
infinities that are caused by integrations over gauge equivalent configurations. 
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1 Introduction 



The standard Faddeev-Popov path integral formulation [0] provides a convenient frame- 
work for calculations in the perturbative regime of QCD. There are, however, some points 
of interest that need clarification: 

• The derivation of the Faddeev-Popov path integral involves infinities from the in- 
tegration over gauge equivalent configurations. Therefore it would be preferable to 
find a derivation of a QCD path integral which is well-defined at all stages. 

• The usual Faddeev-Popov formulation in the continuum gives no prescription for 
how to handle infrared singularities, which may play an essential role for the non- 
perturbative aspects of QCD. 

• There is a lively discussion about the detailed form of the Faddeev-Popov path 
integral, concerning issues like Jacobian factors and the appearance of the Haar 
measure § |, |, |, §. 

In order to clarify these issues we choose canonically quantized QCD as a well-defined 
starting point. By working in 1+1 dimensions we avoid field theoretical infinities indicat- 
ing the need for regularization; only infinities resulting from gauge invariance are left. We 
derive a path integral expression for the transition amplitude in the traditional manner 
by inserting eigenstates of the gauge field operators and coherent states for the fermions 
respectively. With the help of a gauge fixing procedure following Lenz et al. |], we achieve 
formulations purely in terms of physical (i.e. unconstrained) but curvilinear variables. We 
derive two equivalent expressions for the path integral depending on the domain of def- 
inition (compact vs. non-compact) of the remnants of the gauge fields. The inclusion of 
fermions entails no further difficulties and is treated in the appendix. As an illustration 
of the formalism we explicitly calculate the partition function and the spectrum of pure 
SU(iV c ) Yang-Mills theory. Finally we establish the connection to the Faddeev-Popov 
formalism avoiding the usual infinities. 

To make the paper self-contained and to fix our notation we first summarize the formalism. 

2 Review of the Canonical Formalism and Conven- 
tions 

In this section we give a short review of the formulation of 1+1-dimensional Yang-Mills 
theory in the canonical formalism, based on the work of Lenz, Naus and Thies 

The Hamiltonian of 1+1-dimensional Yang-Mills theory in the Weyl gauge is given by 

H= / L dxTr(n(x)) 2 , (1) 
Jo 
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where H(x) is the momentum canonically conjugate to A(x), and equals, up to a sign, the 
chromo-electric field. A(x) and U(x) are su(A r c ) matrices, that can be expanded in terms 
of the N% — 1 standard generators A a /2, 

A(x) = A a (x)y and U(x) = IT(x)y. (2) 

(Repeated indices are summed over throughout the paper.) Space is assumed to be a 
circle, i.e., periodic boundary conditions A(x + L) = A(x) and H(x + L) = H(x) are 
imposed. The theory is quantized by demanding the commutation relations 

[Tl a (x),A\y)} = -i5 ab 5(x-y), (3) 

where 8{x — y) denotes a periodic (^-function. 

In the Weyl gauge, Gaufi's law has to be required as a constraint on physical states, 

Dn|$) = 0, where D ab = 5 ab -^ + gf acb A c . (4) 

Since the operator DII generates small gauge transformations, this means that physical 
wave functions have to be gauge invariant. (In 1+1 dimensions there are no large gauge 
transformations because 7Ti(SU(iV c )) is trivia^].) 

The GauB law constraint is solved in |]J by a transition to curvilinear coordinates (a, A'), 
where the coordinates a constitute a maximal set of gauge invariant quantities, while A' 
represents the gauge variant part of A. Thus, physical wave functions are functions of a 
only. More specifically, a is defined as a constant diagonal matrix containing the phases 
of the eigenvalues of the Wilson loop around the circle: 

e igLa _ y\ V ^ ^ g ^ ^ ^ \ ^ 

where V denotes path ordering and V G SU(A^ C ) is chosen such that 

a = diag(ai, ... ,a Nc ) = a c ° — . (6) 

Because only A^ c — 1 of the iV c variables are independent, a is expanded in terms of the 
N c — 1 generators of the Cartan subalgebra of SU(iV c ) in the last equation^]. 

As can be seen from (|5|), the variables a are defined only up to 

• displacements Oj — > ai H -ki for ki G Z with ki = (7) 



and color permutations Oj — ► o-pa) for P G S 



where S^ c denotes the symmetric group. The set of all displacements will sometimes 
be denoted by V, the set of all color permutations by W (the 'Weyl group' of SU(iV c )). 
Examples for the cases of SU(2) and SU(3) are shown in figs. |l] and ||, resp. 
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Figure 1: Fundamental domain for SU(2). FD denotes the fundamental domain, the 
dashed arrow is a displacement, the dotted arrow (a reflection at the origin) is a color 
permutation 
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Figure 2: Fundamental domain for SU(3). The shaded area is the fundamental domain, 
the dashed arrows are displacements and the dotted arrows (reflections at the diagonals 
of the hexagon) are color permutations 



The invariant hyperplanes of all reflections contained in D x W divide R^ -1 into domains 
in such a way that any a G R^ -1 can be mapped into a given domain by exactly one 
transformation in T> x W. To render the definition of a unique, one can therefore restrict 
its values to one of the domains, which is then called the fundamental domain (FD). A 
special choice for SU(2) and SU(3) is shown in figs. [T] and 

The remaining coordinates A' can be chosen such that the Jacobian of the coordinate 
transformation factorizes, 

= J(a)f[Al (9) 

where J(a) results only from the diagonalization in (|5]) and equals the density of the 

2 Strictly speaking, in the absence of fermions in the fundamental representation the gauge group is 
SU(iV c )/ZAr c , which has a non-trivial fundamental group, but since we are going to include fundamental 
fermions later, we ignore this fact. 

3 The index cq takes on the values i 2 — 1 with i = 2, . . . , N c throughout the paper. 
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'radial' (class) part of the Haar measure of SU(iV c 



In the new coordinates the Hamiltonian ([I]), when acting on physical states, reduces to 
the radial part of the Laplace-Beltrami operator of SU(iV c ), 



H\Q) = H ph \$), 

1 1 d d 

t^t^ — J[a)—— 

2LJ(a)da c ° K ' da c ° 



= — — — • (11) 



We are now in a position to formulate the path integral for this problem. We will derive 
two forms with different domains of a: 

1. a is defined on the covering space R c_1 . 

2. a is restricted to the FD and all integrations are performed over the FD. 

In the first case the transformations ([?]) and (§) will be considered as residual symmetries, 
under which physical wave functions have to be invariant. 



3 Derivation of the Path Integral on the Covering 
Space 



In this section we derive a path integral expression for the time evolution of a physical 
state |<3>), i.e. for the amplitude (A\$(t)), on the covering space. 

As usual, we divide the time interval into a large number N of time slices, each of length 
e = t/N and insert unity at each step, 

r N 

(A\ e- im |$) = / J] [dA-i (A\ e- iH ^|A-i)] (A)|$>, (12) 

i=l 



with the identification An = A. In ( 12 ) H has been replaced by using QTTj ) and the 
fact that H commutes with the projection onto physical states. 

The matrix elements for the infinitesimal time evolution are now evaluated in a 
Schrodinger representation where we pay attention to the fact that the coordinate trans- 
formation yielded a non-trivial Jacobian, 



(A|e _ifr "* e |i4 < _i) = eaq) 



d 2 



IS 



d 



x 



2L \daf> 2 ' JCOK ^da C0 
5(a t - a;_i) 5[A[ - A^_J 



x 



JK) J( Oi _0 JJ'IA^J'IAU 



(13) 



The Jacobians have been split in such a way that we can obtain a simpler form of the 
kinetic energy by moving one of the factors through the differential operators. But before 
this we have to ensure differentiability of the terms the operators act on, which can be 
achieved by use of the identity 



5{a,i 



i-1 



sgn <n sgn Oj_i 



J{a i )J{a i -i) JJ(ai) Jj(ai-x) 



5(di - cij_i), 



(14) 



where sgn a equals +1 if an even number of reflections is needed to map a into the FD 
and —1 otherwise. This sign makes the root of the Jacobian differentiable: 



sgn a J J (a) 



sgn a n 

k<l 



sm 



gL(ai - a k ) 



ocn 

k<l 



sm 



gL(ai - a k ) 



(15) 



If we now move the Jacobian factors to the left, the kinetic energy simplifies and an 
effective potential is induced, 



exp 



d 2 



IS 



2L V <9a c ° 2 



+ /q>(°) 







da c ° 



sgn a 



sgn a 



J(o) Jj(a) 



: exp 
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' 1 d 2 

2Lda c ° 2 



V, 



eff 



(16) 



which is of course reminiscent of the hydrogen atom where V e s would be the centrifugal 
barrier. The procedure just applied corresponds to the transition to "reduced" wave 
functions. Here, however, V e g can be calculated to yield an irrelevant constant 0. It 
follows that 



x 



X 



J(a i )J(a,_ 1 )J'[A'JJ'[^_ 1 ] 
dpi exp 



is \p. 



(17) 



where the integral representation of the 5-function has been inserted and the index cq has 
been supressed from p and a. 

Let us now consider the projection onto physical states. Simply requiring them to be 
independent of the A' with a varying over the whole R 7 ^ -1 is not sufficient to render 
them physical. In addition, we have to postulate invariance under displacements and 
color permutations. A physical state would therefore be, e.g. 

|$) = |a) phys = Yl T a >- (18) 
After insertion into (|12"D this leads to the phase space path integral 4 

. N-l 

/„„,,: , / II da 4 (27r)( 1 - iV ^x 



(Of e 



-iH ph t 



yJj(Of)J(ai) u 



RiV c -l 



i=l 



4 In our final expression for the path integral we perform the projection onto physical states at the 
endpoint. It is irrelevant at which time the projection is done, because the Hamilton operator is itself 
invariant under displacements and permutations. 
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X 



RA/ C 



N 

Yl dpi exp 

i=l 



N 

i£ E (pi 

L 8=1 



tt ;V = Ctf 



(19) 



where we have assumed that <Zf , dj G FD and have performed the integration over the 
unphysical degrees of freedom with help of 



J dA(.) = J daJ{a) J dA'J'[A'}(.). 



(20) 



In eq. ([T9| ) only the sgntu factors of the endpoints are left, while all others have canceled 
pairwise. The configuration space path integral is obtained by performing the Gaussian 
integrations over the momenta, 



(Of e 



-LHphi 



1 



phys 



J(af)J(a,i) 



L s(N c -l)N/2 



2me 



) 



J^sgn 



UJ X 



X 



N-l 



JJ_ daj exp 



i=l 



N 

. i=l 



off 



(2 7\r — tif 



a =a ; 



(21) 



with Of, ai G FD. 



The inclusion of Fermions can be accomplished along similar lines and is therefore post- 
poned to the appendix. 



4 Derivation of the Path Integral on the Fundamen- 
tal Domain 



We start with the Schrodinger equation for the reduced wave function \I/(a) = y J (a) $(a): 

/ ' ^ K ff U(a) =EV(a), (22) 



2L <9a c °< 



together with the condition \I/(a) =0 if J (a) = 0. Because the reduced wave function has 
to vanish on the boundary of the FD, the solution of the Schrodinger equation fl22]) in the 
FD is not influenced by the solutions in its gauge copies. 

We now reformulate the problem. We enforce the vanishing of the wave function on the 
boundary of the FD by introducing a potential, which is equal to zero inside and infinite 
outside the FD. Thus we obtain the new Schrodinger equation 

(~2Zd^ + VcS + Vhox{a) ) * (a) = ^ (G) ' (23) 

with 

T/ . , j : a G FD 

oo : otherwise, 
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and no constraint on the reduced wave function 



As we have a standard kinetic energy and a potential which is bounded from below we 
can apply the Trotter product formula |?J 



lim (e 1 



(V eB +V box (a))t/N 4 -iE kin t/N 



N 



(24) 



section we consider the reduced a-eigenstates J J (a) \a) with the property 



and then derive the path integral in the usual way. To be consistent with the previous 



(25) 



x 



We obtain 



(Of e 



J J(a') (a'\a)J J (a) = 5 (a — a'). 



. N-l . N 

J„_1 »=1 nM.-l »=1 



^Jj(a i )J(a i 
x exp 



RJV C -1 



RiVc-l 



1=1 



2L 



Pi ~ V cS - Vbox(Oi 



Clj\[ — (Xf 



a =a i 



Due to the potential Vb ox , which is rapidly oscillating outside the FD, the a-integration 
can be restricted to give the phase space path integral on the FD 



Of e 



1 



^J(af)J(ai 

x exp 



. N-l . N 

J nda l (2 7 r)( 1 -^^ j ]Jdp. 

FD 1=1 RJVc-i 1=1 

1 



N 



El 0>i O'i—l ' > , - 



& TV — (2f 



(26) 



The configuration space path integral is obtained by performing the Gaussian integrations 
over the momenta, 



(a f | e- lWbox >i) 



1 



^jj{a t )J{a- 1 



.\ Gvrie) 



L ,(N c -l)N/2 



X 



N-l 

> / Y[ ex p 



FD 



1=1 



d]\[=Qjf 

a =ai 



(27) 



5 Equivalence of the Two Formulations 

In the previous sections with have derived two different formulations of the path integral 
for 1+1-dimensional Yang-Mills theory. We will show now that they are equivalent in the 
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limit e — > 0, where we deal with paths which are continuous in time. In the formulation 
on the covering space we then have a sum over endpoints, which are related by small 
gauge transformations and are accompanied by a possible minus sign. We demonstrate 
that only paths which stay within the FD are relevant because all paths that touch the 
boundary or cross it cancel each other. 

To calculate the transition amplitude (of | e _Lffph *|ai) p hy S we have to sum over all paths 
starting at a(0) = a ; and ending at a(t) = w a f including a sum over the symmetry 
transformations uj. 



gLa 




Figure 3: Cancellation of Paths 



Let t\ = inf {r | < r < t, a(r) G <9FD} be the first time that o(r) touches the boundary 
of the FD. Now consider the path 



air) 




0<T<t X 

h<T<t, 



where p G W x T> denotes the reflection at the hyperplane the path a(r) meets at r = t\. 
Because the reflection at the boundary is a symmetry of the Yang-Mills Hamiltonian, 
the two paths a(r) and a(r) have the same action^. But as the final points of the two 
paths differ by a reflection, they have opposite sign and cancel each other in the sum over 
symmetry transformations in eqs. ([19]) and (pip . Thus it is equivalent to consider only 



paths which stay inside the FD and we have arrived at the second form (^6[ [27]) of the 
path integral. 

For the case of SU(2) the cancellation of paths is illustrated in fig. |3|. 

5 This argument holds also for the inclusion of fermions, if all gauge transformations of the a are 
accompanied by the appropriate gauge transformations of the fermions. 
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6 Explicit Calculation of the Transition Amplitude 
and the Partition Function 



To illuminate the role of the sum over endpoints we derive an explicit expression for the 
transition amplitude (pip in terms of energy eigenstates. The Gaussian integrations in 
( 2~I ) can be readily performed to yield the well-known propagator of a free particle plus 
a sum over endpoints, a sign and a Jacobian factor: 



i^^af - a-i) 2 - iV eS t 



(Of e ph ai) P hys = , tt^T 2^ sgnwexp 

yj J (at) J (en) V27rity V 

(28) 

Now we split the sum over u into a sum over color permutations p G W and a sum over 
displacements (JgD. The sum over displacements is turned into a sum over momenta by 
use of Poisson's resummation formula, 

E/(« + *) = ^E [ da'e-^'fia + a'), (29) 



8eV l-^^l p<=X>* 



RiV c -l 



where |EC| is the volume of the elementary cell (the Wigner-Seitz cell of T>) and T>* 
denotes the reciprocal lattice of the Bravais lattice V. It is defined by 



p8 = e 2ttZ for £ G — Z, = 0, (30) 



or, since J^iPi = 0, 



Pi = ^ - -^r E witn n * e Z ' n N c = 0- (31) 

Here, following the convention of (|6]), Pi and 5j denote the diagonal matrix elements of p 
and S, resp. 



Applying (|29|) to fl28 ) and performing the a'-integration we can write 



(a f |e lHph *|ai) phys = E E sgnpexp 

1 J(at)J(a\) |EC| p gd*pgw 



(32) 

The p that are invariant under some permutation p G W, and therefore also under some 
transposition r, do not contribute because for these the exponential is invariant under 
the substitution p — > rp and therefore each term with a given p is canceled by the term 
with rp. For the other p we choose out of each multiplet of N c \ vectors the one with 
Pi > P2 > ■ ■ ■ > Pn c and sum over the permuted ones separately. After rearranging the 
p-sum we obtain 

(«f|e~ i/W lai) phys = 1 E E sgn(pa)x 

/ J{a i )J(ai) | EC | p£V* p.o-ew 

P1>...>PJV C 
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x exp 



(33) 



This formula provides us immediately with the energy spectrum of 1+1- dimensional Yang- 
Mills theory. Using the representation (|31|) and the expression for V e s found in |l[ one 
obtains 



En 



4 



iVc-l 



'JV C -1 



E < - w E 

i=l iVc \ i=l 
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with rii > n 2 > 



and for SU(3) 



> Wjv c -i > 0. For SU(2) this reads 

9 2 L 



E n 



6 



«1 



1) 



(n\ + n\ - nin 2 - 3). 



The corresponding eigenf unctions can also be read off: 

1 



%{a) 



J{a) |EC| p G w 



^ sgnpe ipPa . 



(34) 



(35) 



(36) 



(37) 



These are just the characters of the irreducible representations of SU(iV c ). We see that the 
sum over color permutations together with the sign ensures that the zeros of the Jacobian 
are canceled and the wave functions remain finite. 



Our results agree with those of Rajeev [§], Gupta et al. ||, Engelhardt [|T0[1 , Engelhardt 
and Schreiber [11] and the compact formulation of Hetrick || but disagree with the non- 



compact formulation appearing in the latter. In this formulation the gauge is fixed on 
the classical level and the resulting configuration space, the Cartan subalgebra of SU(iV c ) 
modulo the Weyl group, is quantized under the assumption of a flat measure. As opposed 
to this, in Rajeev's as well as in our (resp. IJ's) approach the full constrained system is 
quantized, and in the course of implementing the constraints the configuration space turns 
out to be the set of classes of SU(iV c ) (the eigenvalues of the Wilson loop) equipped with 
the Haar measure, which is induced by the flat measure of the original (unconstrained) 
configuration space. 

In an earlier work || Hetrick and Hosotani derive the non-compact results also via a 
Faddeev-Popov path integral and find that the Haar measure (the Faddeev-Popov deter- 
minant) is canceled by the v4°-integration. In section |7j we will show that our result is 
equivalent to a Faddeev-Popov-like expression with an additional antisymmetrizing sum 
over gauge copies and that the Haar measure is indeed canceled. The calculation in this 
section shows that these eliminate that part of the spectrum that differs between their 
and our approach. 

The partition function 

Z{0) = Tre -/3H = f da J (a) (a\ e- pH ^ h \a) phys , (38) 

J FT) 
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where the trace is taken only over the physical part of the Hilbert space, contains a matrix 
element that differs from the transition amplitude calculated in section [3] only by a factor 
of i in the exponent. The derivation of (|2T|) and ([33|) is still valid after the replacement 
of it by f3 and \e by e, and so we obtain, after inserting fl3~3f) into fl3~8f) and performing the 
integration, 



E 

pev 
pi>...>pn c 



exp 



2L 



P 2 + V c 



eff 



(39) 



7 Connection with the Faddeev-Popov Formalism 

Having derived the gauge fixed path integral, one is tempted to ask for the connection be- 
tween this approach and the more conventional Faddeev-Popov approach. In this section 
we show that the two expressions for the partition function are actually equivalent if the 
infrared degrees of freedom in the Faddeev-Popov path integral are treated correctly. To 
this end we continue the transition amplitude (^T|) with endpoints identified to Euclidean 
space and insert the result into (j3~8f): 



Z(p) 



r al(o) 

da (0) > sgntu / Va exp 

I D V ^ V MO) 



-/>(^) 2 



off 



(40) 



Here we have adopted the usual continuum notation for the path integral measure and 
integrand. In addition, we have reintroduced the spatial Lorentz index on a. 

The axial gauge corresponds to the gauge condition f[A] = with 



f C0 [A] = <9iA lco + -a 



Oco 



f c [A}=A lc (c^c ), 



(41) 



where a denotes the diagonal part of the spatial zero mode of A . Note that we cannot 
require the zero mode of A lc ° to vanish because of the periodic boundary conditions. 
Instead we use the remaining gauge freedom to eliminate a . 

Now we use ( fill) to introduce the integration over a and the missing modes of A 1 in fl40|) , 



Z(J3) oc J £>a°5>gnw J VA 1 5[f[A}} exp 



A 1 (/3)="A 1 (0) 
o 1 (0)eFD 



p 



dr / dx 



(42) 



where we have suppressed det^i. A 1 contains a 1 as a Cartesian zero mode, and the two 
integrals in (|40f) are contained in the integral over A 1 . In the shorthand notation ^A 1 
the transformation u is understood to act only on the diagonal part of the spatial zero 
mode of A 1 . In the exponent a 1 has been replaced with A 1 due to the gauge condition 
incorporated by the ^-function. 
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As the next step we generate the remaining modes of the time component of the gauge 
field via the identity 



1 = det D; J V'A° exp - jf dr J dx Q (D^ ) 2 + A 1 ^^ 



(43) 



which is valid under the gauge condition. The symbol VA° means that the zero modes 
of the covariant derivative Di are not included in the integral. 

Next we want to show that the Faddeev-Popov determinant is already contained in the 
factors appearing in (|42| ) and ( pE3|) . It can be calculated by looking at the change of the 
the gauge condition under an infinitesimal gauge transformation, 



/3=0 



One obtains 



Sf C0 = -ft (D 1 ^) 00 + -L Jdx (D°5(3) CO , sr = - (D 1 ^)' (c * c ) 



and therefore 



(44) 

(45) 
(46) 



A f [A] oc detD 1 ', 
which is, up to a constant factor, what we found above. 

The full expression therefore reads 

oc^sgnu; JvA"8[f[A\\ A f [A]exp - dr J dx Q (A 1 + D.A 

w A 1 (^)= w A 1 (0) 

(47) 

The above expression differs from the traditional QCD path integral by the sum over 
residual symmetry transformations and the sign factors. It should be noted that these 
features were essential for obtaining the correct spectrum of pure Yang-Mills theory in 
section I5L 



2 _ Kff 

L 



8 Discussion 



In this paper we have derived the path integral for 1+1-dimensional QCD starting from 
canonically quantized QCD in the gauge fixed formulation of Lenz et al. 0. After elimi- 
nating the unphysical degrees of freedom with the help of Gauss's law one is left with a 
set of curvilinear coordinates and the corresponding Jacobian. This yields a non-trivial 
form of the kinetic energy which can be simplified by a procedure analogous to the tran- 
sition to reduced wave functions in the canonical formalism. Depending on the domain of 
definition (compact vs. non-compact) of the remnants of the gauge fields one can derive 
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different expressions for the transition amplitude. We compared the two formulations and 
showed that they are equivalent in the limit of continuous time. We calculated the tran- 
sition amplitude, the partition function and the spectrum in the pure Yang-Mills theory 
and found that the latter agrees with the results of part of the literature f| ||, [|, [Hj |Tl| 
while it disagrees with some other work ||. The connection to the more standard 
Faddeev-Popov formulation could be established by introducing integrations over the 
missing gauge fields. Remainders of the quantum mechanical gauge fixing procedure still 
appearing in the Faddeev-Popov expression are a sum over residual symmetries supplied 
with corresponding sign factors, and an effective potential which is, however, constant in 
this gauge. 

The features which are new to our approach are the direct derivation of the path integral 
from the gauge fixed Hamiltonian instead of the ad hoc introduction of gauge fixing 
terms; the notorious "Faddeev-Popov infinities" are avoided. In our construction residual 
symmetries and boundary conditions on the wave functions are explicitly taken care of. 

The next steps in this direction would be to work out a corresponding path integral 
expression for 3+1 -dimensional QCD. There are problems arising from the fact that the 
path-ordered exponential depends on two space coordinates. 

We hope that the connection between the canonical and path integral approaches provided 
by our calculation can contribute to a better understanding of the renormalization of 
the canonical formulation of QCD. Furthermore, our work might help to clarify the 
relevance of the Jacobian. Because we are using a manifestly gauge invariant approach 
our formulation of QCD is well-suited for all kinds of approximations without the risk 
of violating gauge invariance. One might hope to find traces of non-perturbative physics 
relevant for confinement in 3+1 dimensions in this way. 
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A Inclusion of Fermions 

In this appendix we derive the transition amplitude for the case of fundamental quarks 
coupling to the gauge field. 

The physical Hamiltonian in normal-ordered form (i.e. all ip^ to the left of all ip) reads, 
after the unitary transformation in [|[], 

H ph = H F + H G + H C} (48) 
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where 



and 



Hp = J dxifj^ ^— gaj + 



(3m ] if>, 



(49) 



g 2 f f A a 



-d' 2 



n ab 

{x-y)if> j (x)-^if; l (y) + 



+ 



8 N n 



6N r . 



+ sin" 

k>l 



' gL(ai - a k ) 



(50) 



is the non-Abelian Coulomb interaction. Here, the if) are gauge invariant fermion opera- 
tors, and the form of Hq is that appearing in (|ll|). d' is the covariant derivative in the 
adjoint representation corresponding to a, i.e. d/dx — ig[a, . ], on the space of its non-zero 
modes. The contraction arising in the Coulomb term gives, apart from the term pro- 
portional to the electromagnetic fermionic charge Q m , a term which is proportional to 
the diagonal part of the chromo-electric fermionic charge and vanishes on physical states 
because of the covariant zero mode of Gauss's law: 

Q2l*> = o, (51) 

where |3>) denotes a physical state after the unitary transformation in M. 

We now derive a path integral expression for the amplitude (af,if)f\ e _LH P h *|$), where 
denotes a coherent state for a Dirac spinor. 

The time-slicing procedure is made possible by the identity [O 



-\Ht 



lim 



:e 



AHt/N. 



N 



(52) 



valid for normal-ordered Hamiltonians H. 

In the "time-sliced" transition function we insert at each step the identity in terms of 



coherent states of Dirac spinors [T2| , 

1=1 dijj^dijj exp 



dxif>'if> 



MM, 



(53) 



as well as the identity expressed by eigenstates of gauge- field operators (see fljjD). We 
then evaluate all the operators in a procedure similar to the above. The final result is 



(a N ,iP N \e- lH * bt \$) 

N 

x exp 



t . (JV -l)JV/2 -JV-1 

—J / n [A*M*l>i 



J(ao) sgn q ^ 
\| J{a N ) sgna N 

, ^ . L (ai - ai_i) 2 r M-^i-x, tt tt S 



i=0 



L i=l 

x (a ,ip \$). 



X 



(54) 



In this formula Hp and He depend on the fields if}}, if>i-\ and a;_i. 
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